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INTRODUCTION 

Flutter can result In costly (both In tine and 
money) overruns In turhofan>«ng1ne development pro- 
grams. Solving the problem of flutter (at the 
engine development stage) may well mean major engine 
redesign and retesting. For this reason, engine 
manufacturers and government agencies are currently 
supporting numerous research programs In an attempt 
to develop flutter prediction systems that can be 
used to design flutter-free engines. To date, these 
research programs have Identified five regions of 
the coTiipressor performance map where flutter Is 
general 1y encountered (fig. 1). 

Of the five regions shown In figure 1, the 
supersonic low back-pressure flutter region has been 
the most thoroughly Investigated analytically 
(refs. 1 to 3). In general, these analyses have 
considered the flow field through a cascade of two- 
dimensional airfoils undergoing simple harmonic 
pitching or plunging. The gas stream was assumed to 
be an Inviscid, nonconducting, perfect gas. Shock 
waves that originated In the flow flelo were assumed 
to be weak so that supersonic small -disturbs nee the- 
ory could be used. Although these assumptions ap- 
pear to oversimplify the flow conditions encountered 
by a rotor at the onset of flutter, the flutter 
boundaries predicted by tnese analyses correlate 
well with experimental data. 

A recent analysis (ref. 4) has attempted to 
apply the supersonic, linearized, small -disturbance 
theory to higher backpressure operating conditions 
(I.e., region IV of fig. 1) by Including a finite- 
strength shock wave within the cascade passage. 


Results from this analysis show that the unsteady 
motion of the shock wave tends to Induce bending 
flutter. The existence of this flutter mode Is 
documented In reference 5. 

The remaining operating region where supersonic 
flutter occurs (region V of fig. 1) lies close to 
the stall line of a stage. Analyses of this region 
have not appeared In the open literature. Litera- 
ture on this subject (e.g., refs. 5 to 7) generally 
presents experimental data to document the extent of 
the flutter region. These data provide a very 
limited base from which an empirical correlation can 
be derived for predicting the onset of this flutter 
mode. 

The objective of the present analysis Is to 
develop a model for predicting the onset of super- 
sonic stall flutter encountered by rotors that do 
not have part-span dampers or tip shrouds, kxperl- 
mental data reveal that the flutter mooe is gen- 
erally the first flexural mode at a reduced fre- 
quency (based on tip relative velocity and tip seml- 
choro) of about 0.2. The vibratory pattern around 
the rotor tends to be very regular; All blades 
vibrate at the same frequency but are shifted In 
phase by a positive Interblade phase angle of about 
20 to 50* (ref. 7). This positive phase shift 
Implies that the vibratory pattern Is traveling 
around the wheel In the direction of rotation (I.e., 
a forward-traveling wave). 

A further characteristic of the flow regime 
associated with stall flutter Is shown by the 
steady-state pressure distribution In figure 2. 

This pressure distribution was produced from mea- 
surements taken across a rotor tip while the rotor 
was operating near the flutter boundary. This fig- 
ure clearly Indicates a detached, leading-edge bow 
shock wave Impinging on the suction surface of the 
adjacent blade. The large extent of the compression 
region at the base of the shock wave seems to sug- 
gest a separated flow region that would Increase In 
size as the fan operating point moved toward the 
supersonic stall bending flutter boundary. 
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The present analysis develops a flutter nodal 
for the highly conpllcated flow field Illustrated In 
figure 2 by using tvN>-dlMns1ona1 Ktuator disk the> 
ory. The effects of flow separation are Included In 
the mdel throuoK rotor-toss and deviation-angle 
correlations. For lotr-speed ftowSi actuator disk 
flutter Mdals have bddn able to predict the onset 
of a Slng1e-degree-of-freedo« bending flutter node 
U.g., refs. 8 and 9). The success of the aodets at 
low speeds suggests that a coepresslbte actuator 
disk Model Might be capable of predicting the onset 
of bending flutter at supersonic speeds. This hy- 
pothesis has been conflnaed by a nuMber of calcula- 
tions. These calculations, however, did require as 
Input the Interblade phase angle at the onset of 
flutter. So that this requlrenent could be avoided, 
the actuator disk model was Modified to allow for 
moderate values of Interblade phase angle. This 
modification results In a flutter model that can be 
developed Into a flutter prediction system. The 
validity of this flutter prediction system Is demon- 
strated by comparing the predicted flutter boundary 
of a high-speed fan stage with Its measured boundary. 

FORHULATION 


where the tilde denotes a time-dependent variable 
and the bar slMlfles a steady-state variable. The 
variables a, U(, t, and are the density, 
velocity components, time, and special coordinates. 
respKtlvely. (See fig. a for definition of coordi- 
nate system.) In addition the use of repeated 
Indices denotes summation with respect to the re- 
peated Index. The linearized momentum equations for 
the flow fields are 
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where P*. Is the nondlmenslonal pressure. The 
remaining field equations are the energy equation 
for an Inviscid, nonconducting gas 
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and the equation of state 
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The present analysis for stall bending flutter 
of an Isolated rotor models the rotor as a two- 
dimensional cascade of airfoils. The cascade Is 
defined by the blade-element geometry on a cylindri- 
cal surface at a distance R from the axis of rota- 
tion. The flow field In the cascade plane 1$ 
assumed to be two dimensional, compressible, and 
time dependent. Viscous forces are considered only 
within the blade channel. The unsteady flow vari- 
ables associated with the rotor vibratory motion are 
assumed to be smaller than their steady-state coun- 
terparts. These variables at an Instant In time are 
required to be periodic around the wheel at a period 
equal to a fractional part of the circumferential 
distance d ■ 2iR. 

The motion of the aWolls In the cascade plane 
Is restricted to simple harmonic plunging and edge- 
wise motion at a cyclic frequency « (fig. 3). In 
addition the motion of each airfoil at an Instant In 
time Is assumed to be shifted from that of Its 
neighbor by an interblade phase angle of o ■ 2<n/N, 
where n Is an Integer and N Is the number of 
blades In the rotor. In the present analysis the 
reduced-frequercy parameter kk and the Inter- 
blade phase angle are assumeo to be small. (Reduced 
frequency k(, - where b Is the semi- 

chord of an airfoil and '5_„ Is the magnitude of 
the relative Inlet velocity In the cascade plane.) 

Field Equations 

I he governing linearized equations for the flow 
field upstream and downstream of the cascade are 
written with respect to the relative coordinate sys- 
tem. The variables In these equations are wn- 
dlmenslonallzed by the steady Inlet speed 
the circumferential distance around the wheel d, 
the Inlet steady-state fluid density ana the 
Inlet steady-state fluid temperature T^. In 
addition the subscripts (>, —) signify downstream 
and upstream flow variables. 

The linearized continuity equation for the up- 
stream and downstream flow fields Is 
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The variables appearing In equations (3) and ( 4 ) are 
the entropy s; the ratio of specific heats v; and 
the Inlet steady-state relative Nach number ICm. 
(Theentropy was nondlmenslonal 1 zed with respect to 

t_/«^) 

The present analysis assumes that there are no 
sources of entropy or vorticity upstream of the cas- 
cade. The upstream, unsteady velocity field must 
therefore be Irrotatlonal and hence Is equal to the 
gradient of a potential function a^, where 


“ 1 ,- - 1 - 1. 2 (5) 

The unsteady velocity field downstream of the cas- 
cade Is expressed as 


“ 1 .- " Uq ^ 1-1.2 (6) 

where aiJaXi represents the Irrotatlonal compo- 
nent of the field and H . Its rotational 
component. The source of the rotation field Is the 
vorticity shed by the oscillating airfoils. 

The solution to the field equations Is obtained 
by assuming a simple harmonic spaclal and temporlal 
dependency for the potential, the rotational veloc- 
ity fields and the downstream entropy field. The 
details of the solution procedure can be found In 
reference 10. 

Boundary Conditions 

It Is shown In reference 10 that four boundary 
conditions must be specified at the actuator disk to 
relate the unsteady flow field to the cascade defor- 
mation and the steady Inlet flow conditions. 

The first of these boundary condition requires 
the flow to be continuous across the deforming 
disk. The analytical form of this boundary condi- 
tion Is 
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where represents the nomal velocity 

components relative to the defonnlng disk. These 
velocity components to first order (I.e. . small 
cascade-plane deflection) are equal to 


a - s' 

•• Hlj 

The variable In this equation Is the angle 
between the chora line and the tangent to the camber 
line at the leading edge. The corresponding rela- 
tion for the local deviation angle Is 
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where h^ Is the displacement of the dlsr from 
Its mean position In the axial (Xi) direction 
(fig. 4). 

The next boundary condition requires the total 
enthalpy with respect to the defonnlng actuator disk 
to be locally conserved across the disk. The equiv- 
alent mathematical statement Is 
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where el Is the 1oca1» relative exit flow angle 
measured from the normal to the tral ling-edge plane 
and Is the angle between the chord line 
and thenangent to the camber line at the trailing 
edge. 

Relationships for the angles. Inlet Mach number 
and cascade solidity in the above equations are 
derived by expanding these variables along with 
equations (13), (14) and (lb) to first order. The 
details of this expansion can be fiund in refer- 
ence 10. 

The last boundary condition specifies the local 
entropy tnat Is generated as the fluid passes 
through the cascade. The loss associated with the 
local entropy rise Is defined In terms of a loss 
coefficient defined as 


T-Kjljr T. atXj.O (10) 

where T«. Is the upstream and downstream gas 
temperature and qSS? is the magnitude of the 
relative velocity approacring and leaving the disk. 
The expressions for qlv are 
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The variable in this 

equation represents the local total-pressure loss 
measured relative to the deforming actuator disk. 
The entropy rise across the disk is related to this 
quantity by the equation of state 



Thus s. must be equal to 


(17) 



where h? is the local displacement of the cas- 
cade in the tangential (Xgj direction (fig. 4). 

The third boundary condition requires that the 
local inlet and exit flow angles relative to the 
deforming disk satisfy a deviation angle correlation 
of the form 

a « g(a, u, e‘ ) ( 13) 


The analysis assumes that the local loss co- 
efficient X is related to the local inlet flow 
field, as observed with respect to the deforming 
disk, by a relationship of the form 

X ■ x(a, M__, M, e' ) ( 18) 

As in equation (13) the parameters appearing in this 
equation can be related to the unsteady flow fi>1d 
approaching the cascade and the cascade aefrT<i- 
tion. Upon combining equation (17) with equation 
(18) the entropy rise across the cascade can be 
evaluated in terms of these perturbed quaniiles. 


where o is the local, relative incidence angle, 
is the local, relative inlet Mach number, u 
Is the local solidity of the deformed cascade, and 
e' is the local stagger angle measured from the 
normal to the leading-edge plane (fig. 4). The 
local, relative incidence angle is related to the 
local, relative inlet flow angle el. measured 
from the normal to the leading-edge plane by the 
equation 


Aerodynamic Force 

Having described the procedure for ottermining 
the flow field surrounding the cascade, the unsteady 
forces acting on the cascade and hence its stability 
can be assessed. The aerodynamic force acting on 
the cascade is obtained by considering the flow 
field through ^ control volume that is fixed to a 
cascade passage (fig. 3). The present analysis 
assumes that the relative motion between neighboring 
airfoils, a measure of which is the interblade phase 
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•ngl«, 1$ Mtll and that the reduced frequency of 
thu notion 1$ also snail. The first assunptlon 
Inpllas that the spatial variations of the flow 
variables across a cascade passage are snail and 
thus can be neglected. The low-frequency assunptlon 
suggests that the rate of change of nass and nonen- 
tum within the control volume Is smaller than tne 
net mass and momentum flux across the control volume 
surface and thus will also be neglected. Based on 
these considerations the llnearlied nomentum equa- 
tion for the control volume illustrated In flgur. 3 
can oe written as 


• (P, - P.,)~ ♦ (P« - P.J^ 


♦ - ITj ♦ ff(Uj 
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where F, Is tne force exerted by the airfoils on 
the control volume, t Is the pitch of the cascade, 
U**®' IS the mass flow exiting the control volume, 
and Is the unit vector normal to the exit 
plane of the control volume. (This equation can be 
derived by expanding the quasi-steady momentum equa- 
tions for the control volume shown In figure 3 to 
first order in perturbed variables. The details can 
be found In reference IG.) 

The work done by the aerodynamic force on an 
airfoil over a cycle of motion Is equal to 



where fig [ ] denotes the real part of the expres- 
sion and the asterisk superscript signifies the com- 
plex conjugate of tne variable. The airfolTs 
motion can be described by 
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where the damping coefficient C Is defined as 


C • ■jjS k (24) 

The variable a In this equation Is the log decre> 
ment of the mechanical damping, and N5 Is the 
mass per unit length of a rotor blade. This vari- 
able Is defined as 


• Ko„tc (25) 

where K Is a constant of proportionality, »|| Is 
the blade density, t Is the average thickness of 
the blade, and c Is Its average chord length. 
Introducing this equation along with equation (24) 
Into the inteir jted form of equation (23) yields 

Equating equation (26) with equation (22) and solv- 
Ing for a establishes the minimum level of mech- 
anical damping required for stability at an operat- 
ing point as 

(2T) 


The symbol Mlni^. denotes the mlniiMim value of 
the function with respect to kp. If the available 
mechanical damping of a rotor blade exceeds this 
critical level, any small bending motion Imparted to 
the blade will decay In time. Hence the system Is 
stable and flutter will not occur. However, If 
^Onsct greater than the available mechanical 
damping of a rotor blade, aqy small binding motion 
Imparted to the blade will cause the blade to flutter. 

RESULTS 


where k^, the wave number in the Xj direc- 
tion, 1» equal to k^ - 2m. Equation (21) can be 
Integrated to yield 

“Aero • iN) 

where L J represents the imaginary part of the 
expression within the brackets. At the onset of 
flutter the aerodynamic work per cycle Is equal to 
the mechanical energy dissipated over a cycle, if 
the mechanical dissipative force Is proportional to 
the velocity of the airfoil, the mechanical energy 
dissipated over a cycle of motion Is 


For the flutter mode under study the aaqilltude 
of the motion increases monotonically along the blade 
span from tne node line. For blades that are rigidly 
fixed at their root (I.e., mechanically constrained), 
the node line can be assumed to lie outboard of the 
blade platform (I.e., outboard of the aerodynamic hub 
of the blade). The node line of the first flexural 
mode of the rotor to be investigated later In this 
section Is approximately 20 percent of span height 
Outboard of the platfcrm. Thus, the vibratory motion 
of the blade in the rotor at any spanwise section Is 

hj . -t sin (28) 


h^ - t cos 129) 
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Mhert { it ^ Mpittudt of th« aotlon at tiM sec- 
tion end «| the stagger angle at 20 per- 

cent of span. 

The aerodynamic work per unit of span at a 
given radial location depends on the «g>11tude of 
the motion of the section and the relative dynamic 
pressure of the Incoming streamline to the section. 
Thus hecause both Increase with dlstaire from the 
hub. the total work done by the airstream on a blade 
Is strongly Influenced by the unsteady flow field 
surrounding the tip region. A simple calculation 
shows that this Influence Is concentrated In the 
outer 2S percent of the span. If the outer-casing 
boundary layer Is assumed to Influence the flow 
field over S percent of the tip region of a blade, a 
direct correlation might exist between the flutter 
boundary of a rotor and a twt-dlmenslonal cascade 
whose geometry and dynamic response coincide with 
that of the rotor at 85 percent of span (I.e., alge- 
braic mean of 75 and 95 percent of span). The re- 
sults presented In the remainder of this section are 
based on this premise. 

The objective of the first series of calcula- 
tions Is to establish the Influence of reduced fre- 
quency kk and Interblade phase angle on the 
aeroelastic stability of a cascade of airfoils. For 
this study a normalized damping parameter defined as 



Is calculated as a function of k|, and the Inter- 
blade phase angle. The parameter cq in this 
equation Is the Inlet stagnation density measured In 
the absolute coordinate system. For in greater 
than zero the airstream Is supplying energy to the 
cascade. If the mechanical damping of the cascade 
system Is zero, the operating condition would be 
unstable. For values of less than zero the 
cascade Is doing work on the airstream, and hence 
the system Is stable. 

Mlving equation (30) requires as Input Infor- 
mation the steady-state Inlet flow properties, the 
geometry of the cascade, and the frequency of os- 
cillation. A set of loss and deviation-angle cor- 
relations must also be specified. The present 
analysis makes use of the correlations derived In 
reference 10. The steady-state Inlet flow condi- 
tions are derived from the blate-element data, mea- 
sured at 65 percent of design speed, of the second 
stage of the NASA 14SO-ft/sec-t Ip-speed two-stage 
fan. This blade-element data set Is reported in 
reference 5. The cascade geometry represents the 
geometry of the second-stage rotor element at 85 
percent of span height from the hub. 

Figure 5 shows a plot of aw * function 
of Interblade phase angle and reouceo frequency 
kK for the second stage rotor operating at the 
flutter-boundary at 85 percent of design speed. 

This figure shows that the bending flutter mode Is 
associated with a positive Interblade phase angle 
(I.e., Implies that the vibration pattern Is travel- 
ing around the rotor In the direction of rotation). 
The figure also shows that Increasing the reduced 
frequency stabilizes the motion at a constant tnter- 
blade phase angle. 


Figure 6 shows a plot of «s a function 
of Interblade phase angle and percentage of design 
weight flow along the 65-percent-speed line. For a 
olven Interblaoe phase angle, decreasing the weight 
flow at constant wheel speed (I.e., Increasing the 
steady-state aerodynamic loading) tends to 
destabilize the rotor. 

The results piesented In figures 5 and b are 
also consistent with the observed characteristics of 
supersonic bending flutter reported In reference 7. 

A deficiency of the present theory, as shown by 
the results In figures 5 and b. Is that It predicts 
the damping, to Increase monotlonically with 
Interblade phase angle. This trend Is physically 
unrealistic, and Is a consequence of the small 
Interblaoe phase angle assumption Inherent In actua- 
tor disk theory. To overcome this deficiency an 
additional assumption was Incorporated Into the the- 
ory which limited the value of the aerodynamic damp- 
ing. The assumption is made to replace o by sin 
a everywhere In the analysis. This assumption Is 
mathematically constant with the actuator disk 
model. Figure 7 shows the results obtained with 
this assumption Incorporated Into the analysis. 

For |o| < 45*, the modified analysis yields results 
that are comparable to those obtained from the 
original formulation. In addition the modified 
analysis predicts a finite maximum value for aq 
for 0 • 90*. This maximum value for damping Is 
equal to the value computed at a • 1 radian using 
the original formulation. Experimental results of 
reference 7 show that this value for a (I.e., a ■ 1 
radian) corresponds rather closely to the measured 
Interblade phase angle of the least stable flutter 
mode. 

An additional validation of the modified analy- 
sis as a flutter prediction system for supersonic 
stall bending flutter can be established by showing 
that, for a given rotor, flutter will oe observed 
whenever the predicted xiaximum value of am ex- 
ceeds the normalized structural damping of a rotor 
assembly. For this demonstration the measured flut- 
ter boundary of a scaled model of Fan C of the NASA 
Quiet Engine Program (ref. 6) Is correlated with the 
predicted boundary. The performance map for the 
scale model (tig. 6) shows two flutter Instability 
zones. The first zone is torsional flutter that 
occurs at part speed and lies close to the stall 
line. Just above this zone lies a zone of bending 
flutter that appears to extend past 95 percent of 
design speed. Only the bending flutter mode of this 
rotor is analyzed in this paper. 

Figure 9 shows a plot of the phase shift be- 
tween the amplitude time history of a rotor blade 
and the reference blade No. 17 at the 95 percent of 
design speed flutter point. The data Implies the 
flutter mode is predomlnatly a 3 nodal diameter for- 
ward traveling wave which corresponds to an inter- 
blade phase angle of approximately 41*. This value 
for the interblaoe phase angle of the flutter mode 
correlates well with the data of reference 7 and 
lies within tne unstable range predicted by the 
analysis. Figure li) shows the calculated maximum 
normalized aerodynamic damping for this rotor as a 
function of weight flow and wheel speed. The re- 
sults presented on this curve were calculated using 
the blade element data recorded at tne five operat- 
ing condition shown on figure 8. TV shoded region 
on figure 10 represents the unstaole t utter 
region. The boundary of tnis region ..i'. estimated 
by assuming a log decrement of O.U<^5 fur the mech- 
anical damping of the assembly (ref. 3). Hence all 
operating conditions whicn have an aerodynamic damp- 
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Ing level In excess of this value are assuMd to be 
within the flutter region, while those that have an 
aerodmamlc damping level less than 0.02!t are stable. 

The results on figure 10 clearly show the 
significant Influence of wheel speeo and aerodynamic 
loading on the aerodynamic damping. The calculated 
aerodynamic damping at 60 and 70 percent wheel speed 
lie significantly below the Instability boundary. 

No bending flutter was observed at tnese operating 
conditions. At 90 percent of deslon speed bending 
flutter was observed near the stall region. The 
calculated aerodynamic damping for the near stall 
operating point Is seen to He just Inside the un- 
stable region. Bending flutter was also observed at 
95 percent of design speed near the stall boundary. 
The calculated aerodynamic damping for this operat- 
ing point lies well within the ui>stable region of 
figure 10. At the Internedlate operating point on 
the 9s percent speeO line the calculativ* aerodynamic 
damping falls below tllr Instability bour^lary. 
Therefore the predlcteo flutter boundary at this 
wheel speed lies between the measured flutter bound- 
ary and the Intermediate operating point. Tran- 
scribing the stability boundary of Hgure 10 onto 
the Fan C performance map yields a theoretical flut- 
ter boundary, which Is shown as a dashed line on 
figure 11. The measured operating conditions at 
which bending flutter occurred are shown on this 
figure as solid symbols. The overall agreement 
between theory and measurement Is very good, with 
the analysis slightly over predicting the extent of 
the flutter boundary at 95 percent of wheel speed. 
Below 90 percent of wheel speed the analysis pre- 
dicts the flutter boundary will bend back Into the 
stall regiun. This result is confirmed by the 
experimental measurements. 

The results presented in this section clearly 
show that the current analysis can oe used to pre- 
dict the onset of supersonic bending flutter In un- 
shrouded fans. The input variables to the analysis 
can all be derived from fan structural and aero- 
dynamic design variables. The governing equations 
are simple algebraic equations, and hence are easily 
programmed. These features provide the designer 
with a simple and reliable model for analyzing the 
susceptibility of a fan design to supersonic stall 
bending flutter. In addition it provides him with a 
tool that can also be used to evaluate proposeo 
fixes to an existing flutter problem. 

The engineering approximation introduced to 
extended the range of validity of actuator disk 
theory appears reasonable. To develop a correct 
first principles extension of actuator disk theory 
to finite interblade phase angle for the flow 
conditions as they exist at the tip of a rotor at 
the onset of supersonic stall flutter would re- 
quire a viscous transonic cascade analysis for 
finite interblade phase angle. No such analysis 
appears tu be forth coming in tne near future. Any 
extension which is derived from inviscid aerodynamic 
theory would be an approximation, tne result perhaps 
being no more accurate then the simple extension 
developed in this analysis. 

CONCLUDING REMARKS 

In the present investigation an analytical 
model was developed to predict the onset of siver- 
sonic stall bending flutter in high-speeo rotors. 

The stability boundary predicted by the analysis was 
shown to be in good agreement with the measured 
boundary of a high speed fan. The analysis also 
predicted that the flutter mode would be a forward 


traveling wave, be sensitive to wheel speed and 
aenxtynamlc loading. This was confirmed hy the 
experimental measurements. In addition the analysis 
showed that reduced frequency and dynamic head also 
pl«y a significant role In establishing the super- 
sonic stall bending flutter boundary of an un- 
shrouded frn. 

This analysis provides the designer with a sla»- 
ple and efficient prediction tool which he can use 
to analyze tne susceptibility of a fan design tn 
supersonic stall bending flutter. 

APPENDIX A 

SYHBCX.S 

b airfoil semichord 

c average blade chord and damping coefficient 

d circumferential distance around wheel 

F j aerodynamic force 

h^ axial dlsplKoment of an airfoil 

Ii2 tangential displacement of an airfoil 

1 sA 

K constant 

k reduced frequency ba&ed on circimiferentlal 

distance, wd/q^ 

k|, reduced frequency based on semichord, 

wb/q_^ 

wave number of disturbance along cascade 
plane 

M Nach number 

mass of a rotor blade per unit length 
N number of blades In rotor 

n Integer 

n^ unit vector 

P total pressure 

p static pressure 

q magnitude of velocity field 

R radial distance 

s entropy 

T temperature 

t time 

velocity components 

relative velocity component normal to disk 
rotational velocity 
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Figure 1. - Compressor performance and stability map. 
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gure 2. - Typical tip rotor static-pressure distribution near supersonic stall 
flutter boundary. 
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Figure 4. - Geometry of deformed cascade. 


Enlarged view 

Figure 3. - Cascade plane deformation. 
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Figure 5. - Calculated, normalized damping as a (unction of Interblade 
phase angle and reduced frequency for NASA two-stage (aa 
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Figure <c - Calculated, normalized damping as a function of 
interblade phase angle and weigM flow for NASA (wo-stage fan 
at 8S percent of de^n edteel spaed. 
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Figure 7. - Original and modified analyses for NASA heo-stage 
fan at 85 percent of design wheel sp^ and 75. 88 percent of design 
weight flow. 
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Figure 1 - NASA ntoOng rigt tei C model 
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Figure 10. - Maximum normalized damping as a function of weight 
flow and wheel speed for fan C model. 
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figure IL - NASA rotating rig, fan C model. 


